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ABSTRACT 

Mathematical models are developed for the hydrodynamic 
and aerodynamic forces which act on a sailing boat. The 
solution of the force and moment equilibrium provides the 
sailing boat's speed and attitude for performance prediction. 
Each force model is evaluated for its capability to make 
adequate engineering predictions. Further, expressions are 
given for the dimensions and coefficients to provide a feasible 
baseline design to which parametric variations may be made. 

The system of force equations constitute a parameterized 
system model in four degrees of freedom for steady state 
sailing conditions. The mathematical model provides a design 
tool for the naval architect conducting geometric parameter 
trade-off analysis against hydrodynamic performance. 

The principal result is a computerized model which is 
capable of predicting boat speed, heel angle, leeway angle 
and rudder angle for a geometrically defined sailing boat in 
a given wind condition. 'Example boats are evaluated and the 
results recorded in order to demonstrate the use and versa- 
tility of the model. 
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NOMENCLATURE 



A k 

A 1 

A m 



A 

r 

A 

s 

A t - 

BWL 

C dc = 
C do = 
C dis - 
C f = 

C ffb 
C fk 

C fm = 

C, 

fr 

C. 

lr 

C. 

Is 

C mc/4 - 
C 

n 



C popt 



keel cross-sectional area. 

projected lateral area of the hull above the LWL. 

cross-sectional area of the maximum hull station, 
or cross-sectional area of mast depending on 
context. 

platform area of the rudder, 
total area of the sails. 

projected transverse area of the hull above the LWL. 

maximum waterline beam. 

rudder cross flow drag coefficient. 

rudder section form drag coefficient. 

sail induced drag coefficient. 

friction drag coefficient. 

fairbody friction drag coefficient. 

keel friction drag coefficient. 

hull form drag coefficient. 

rudder friction drag coefficient. 

rudder lift coefficient. 

sail lift coefficient. 

rudder moment coefficient about the mean quarter 
chord. 

rudder normal force coefficient, 
prismatic coefficient, 
parent prismatic coefficient, 
optimum prismatic coefficient, 
rudder drag coefficient. 
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c 

y 


keel side force coefficient 


C ysb ~ 


slender-body side force coefficient. 


n 

n 

ii 


rudder side force coefficient. 


D s = 


sail drag 


D so = 


offwind sail drag. 


F h = 


aerodynamic hull force. 


F 

n 


Froude Number. 


H 


maximum canoe body draft. 


K 


heeling moment. 


K h = 


hydrostatic righting moment 


K k = 


keel heeling moment. 


K r = 


rudder heeling moment. 


K s = 


sail heeling moment. 


L 


load waterline length. 


L s = 


sail lift. 


L so = 


offwind sail lift. 


LCB 


position of center of buoyancy (%LWL) . 


LCB 

o 


parent position of center of buoyancy. 


LCB = 

opt 


optimum location of center of buoyancy. 


LWL 


load waterline length. 


M 


pitching moment 


N 


yawing moment 


N k = 


keel yawing moment. 


N 

r 


rudder yawing moment. 


N sb = 


slender-body yawing moment. 


N sx = 


yawing moment due to sail driving force 


N sy = 


yawing moment due to sail side force. 
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p 



vector of assumed variable values. 



P l,2,3, or 
R f = 

R fm = 



R 

w 

R 

wp 

R 

R-, ^ ^ 
1,2,3, or 

s fb = 

s k = 

S r = 

T c = 

u b = 

V a = 

V b = 

V t = 

X 

X h " 

X k = 

X 

o 



Y k = 



^assumed value of individual variable. 

total frictional resistance. 

total form resistance. 

Reynold's Number. 

upright resistance. 

wave resistance. 

parent wave resistance. 

vector of force equilibrium errors. 

individual force equilibrium errors. 

fairbody wetted surface. 

keel wetted surface. 

rudder wetted surface. 

canoe body maximum draft. 

boat velocity in the X q direction. 

apparent wind velocity. 

boat velocity parellel to the y axis. 

true wind velocity. 

force in the X direction, 
o 

hull resistance. 

keel and hull induced drag. 

axis on the plane of the free surface in the 
direction of boat velocity. 

rudder induced drag. 

sail driving force. 

force in the Y q direction, side force, 
keel side force. 
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Y 

o 


an axis on the plane of the free surface at 
right angles to boat velocity. 


Y 

r 


rudder side force. 


Y 

s 


sail side force. 


Z 


force in the Z q direction, heave force. 


Z 

o 


an axis perpendicular to the free surface 
downward . 


a f = 


fortriangle aspect ratio. 


a m 


mainsail aspect ratio. 


a r 


rudder aspect ratio. 


a = 

s 


total sail area aspect ratio. 


b b 


draft at forward perpendicular. 


ii 

0 

XI 


maximum draft. 


b r = 


rudder span. 


b l,2,3 = 


intermediate draft measurements. 


c rr = 


rudder root chord. 


c rt 


rudder tip chord. 


II 

u 

u 


rudder quarter chord. 


O 1 
II 


keel mean chord. 


c = 

r 


rudder mean chord. 


n 

• — i 


metacentric height/LWL. 


d 2 “ 


non-linear hydrostatic righting moment coefficient. 


d 3 = 


Froude Number dependent righting moment coefficient 


II 

•H 


spline cubic function. 


g 


acceleration due to gravity. 


r o 


canoe body draft at maximum draft measurement. 


s = 


sinkage . 


s ij = 


9R./3P • • 
1 J 
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ft 

* 

II 


taper ratio. 


X = 


axis in the load waterline plane along the 
centerline of the boat. 


II 

rQ 

X 


horizontal distance to the maximum draft 
measurement. 


x cph 


horizontal distance to the aerodynamic center 
of effort on the hull. 


x cpk 


horizontal distance to center of effort on 
the keel. 


x cpr 


horizontal distance to the center of effort on 
the rudder. 


x cps 


horizontal distance to center of effort on 
the sail. 


x e 


main boom length. 


II 

•m 

X 


base of the fortriangle. 


x m 


horizontal distance to the mast. 


x = 

o 


0.0, location of maximum draft. 


x pc/4 


horizontal distance from the rudder mean chord 
quarter point to the rudder center of effort. 


X = 

rs 


horizontal distance to the rudder stock. 


X 1 , 2 , 3 = 


horizontal distance from maximum draft to 
intermediate draft measurements. 


y 


axis orthaginal to the x axis to starboard 
in the waterline plane. 


Z = 


axis orthaginal to the waterline plane, downward 


z cpk - 


vertical distance to the center of effort on 
the keel. 


z cpr 


vertical distance to the center of effort on 
the rudder. 


z cpxs 


vertical distance to the sail driving force. 


z cpys 


vertical distance to the sail side force. 


z fb 


freeboard. 



12 




z . 
1 



z 

p 



z 



ro 



z = 

rrc 

V 

c 



A 



a — 

r 

6 

e u = 

Y 

6 

0 

X 



Y h 2 o 

TT 



P • 

K air 

P h 2 o - 

<p 

ip 






vertical distance to the center of gravity, 
mast height, 
main hoist. 

vertical distance from the rudder root chord to 
the center of effort on the rudder. 

vertical distance to the rudder root chord. 

volumetric displacement of the fairbody. 

volumetric displacement of the keel. 

weight of the boat in tons. 

rudder angle. 

apparent wind angle. 

upright apparent wind angle. 

true wind angle. 

variable increment. 

pitch angle. 

angle between boat centerline and aerodynamic 
hull force. 

kinematic viscosity of water. 

3.14159 

density of air. 
density of water, 
heel angle, 
leeway angle, 
rudder sweep angle. 
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I. INTRODUCTION 



1 . 1 General 

All designers of engineering systems desire a procedure 
for predicting the performance of their designs prior to 
final construction. Until lately, prediction of a sailing 
yacht's performance has been limited to evaluation of scale 
models and intuition. The development of the high speed 
computer and numerical techniques for rapid solution of non- 
linear equations has made the use of mathematical models a 
reasonable alternative for performance prediction. It is 
the purpose of this paper to present mathematical models for 
the forces and moments acting on a sailing boat. The solution 
of the force and moment equilibrium will provide the boat's 
speed and attitude for performance prediction. 

Each force is parameterized in dimensions and coefficients 
generally utilized by designers, to make the model a useful 
design tool for the naval architect. Further, expressions 
are given for the required dimensions and coefficients to 
provide a feasible baseline design from which parametric 
variations may be made. Each force model is evaluated for 
its capability to make adequate engineering predictions, 
including inherent limitations. 

1 . 2 Historical Development 

From the earliest days, yacht design has been called an 
"art". Though the use of wind propulsion for ships is ancient 
the action of the forces and moments on the sailing ship is 



complex and multidimensional. Historically there has been 
no technique for simple evaluation of the effect a change in 
one parameter has on the forces acting on the boat. 

Until fifty years ago design evolved slowly through ober- 
vation of existing ships. Due to the inability to evaluate 
a change prior to incorporation into a newly constructed 
yacht, new designs tended to consist of only conservative 
changes to a previously successful design. This was under- 
standable for failure of an innovation in the final yacht 
could be professionally catastrophic for the designer. 

In the early 1900 's, techniques were developed for pre- 
dicting yacht performance from model tests. A designer could 
now build a relatively inexpensive model and evaluate it 
against other model results prior to incorporating a modifica- 
tion into a new design. Adequate prediction of full scale 
results have been limited by adverse scale effects, but the 
usefulness of these tests is unquestioned. Unfortunately, 
the cost of model testing still precludes extensive parametric 
trade-offs. Each parametric change requires construction and 
testing of individual models. Usually, the limited design 
budget for one yacht precludes any systematic variation of 
parameters, and testing is limited to one or two candidate 
models. The appeal of a parameterized mathematical model 
which is inexpensive to utilize is apparent. 



15 



The idea of using mathematical expressions for describing 
yacht performance is by no means a new concept. The rating 
rules used to handicap yachts are an attempt to apply 
mathematical models to boat performance. Because rating 
rules attempt to measure a boat's speed potential for a large 
range of conditions and reduce the results to a single number, 
it is easily understood why they have little application to 
predicting performance for a given wind direction and speed. 

It was not until the second half of this century that 
towing tank test results were combined with sail force 
coefficients in a computer to solve force and moment equili- 
briums. These programs could provide prediction of performance 
for a given condition within the limitations of the data. 

These solutions were for the specific boat tested and a 
scaled family of geometrically similar boats. 

In the early 1970 's attempts for steady state solution 
from mathematical models began to be published. One pioneer 
in the field has been John S. Letcher, publishing parts of 
his model in references (1) , (2) , and (3) . Further work has 

been done by Myers (4) and Kerwin (5) . The computational 
technique for solution of the force equilibrium condition 
has been described by Hazen (6) . 

By parameterizing forces in four degrees of freedom, the 
model developed in this paper grows from these previous works. 
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1.3 Overview 



In order to resolve the aerodynamic, hydrodynamic and 
hydrostatic forces acting on a sailing boat a coordinate 
system and origin are established. Achieving an equilibrium 
of forces and moments requires that the summation of all 
forces and moments acting either through the origin and about 
the origin are equal to zero. To facilitate modeling, 
component forces which act independently or nearly so are 
identified and examined. Once the forces acting on the system 
and the location and direction through which they act are 
identified, each may be modeled and superimposed to give 
the required equilibrium solution. 

1.3.1 Coordinate System 

Figure 1 presents a coordinate system which describes 

the sailing boats velocity and attitude. The coordinate 

system designated by the capital letters X Q , Y , and Z q are 

located on the free surface with X and Y axis in the 

o o 

horizontal plane. The X q axis lies in the direction of the 

boat's motion and the Y axis is to starboard at right angles 

o 

to this velocity. The Z q axis is vertically downward. 

A second coordinate system is fixed in the boat at the 
forward end of the waterline as presented in Figure 2. This 
coordinate system is designated by the small letters x, y, 
and z. The x axis is oriented along the fore and aft center- 
line on the boat's waterline with the y axis in the transverse 
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Figure I - Coordinate Systems 
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Figure 2 



Boat Coordinate System 
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direction to starboard. The z axis desends orthogonally to 
the x-y plane in the direction of the keel. 

The angular differences between the two coordinate 
systems describe the attitude of the boat as it moves through 
the water. The positive direction for these angles are also 
presented in Figure 2. Their direction was selected to pro- 
vide positive values when the boat is in equilibrium under 
normal sailing conditions with a weather helm. The order in 
which these angles are taken is important to the final atti- 
tude. The boat is first rotated about the Z q axis through 
the angle ip . The angle ip is identical to the boat's leeway 
angle. The boat is then rotated about the X axis, an angle <p , 
which is the boat's heel angle. The boat is then rotated 
about the Y q axis, an angle 0 . 9 is the angle by 

which the boat pitches down at the bow and is also the angle 
of the axis about which the boat heels. Finally the angle a R 
is the angle the rudder is deflected from the boat's centerline. 

A sailing boat responding to the forces acting upon it 
and conforming hydrostatically to deformations in the free 
surface would heave in the positive Z q direction an amount s. 

The selection of these coordinate systems allows the 
modeling of forces in either coordinate system. When modeling 
velocity related forces, such as wave drag and lift, the 
coordinate system fixed in the water is more attractive 
computationally. On the other hand, hydrostatic forces are 
more easily computed in the coordinate system fixed in the 
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boat. In the final solution of force equilibrium, forces 
and moments are resolved in the X , Y q , and Z q system except 
the heeling moment which is resolved about the x axis. 

It is a common practice in naval architecture to describe 
a vessel's dimensions and coefficients in a positive sense. 
With few exceptions this is done throughout this paper. 

The positive direction of forces and moments is defined in 
Section 1.3.2. 

1.3.2 Equilibrium of Forces 

The forces and moments acting upon a sailing boat 

in these coordinate systems may be divided into six components 

X = Forces in the positive X Q direction 

Y = Forces in the positive Y q direction 

Z = Forces in the positive Z q direction 

K = Moments about the x axis 

M = Moments about the y axis 

N = Moments about the z axis 

o 

The positive direction of the moments K, M, and N are taken 
in the same sense as the angles <f> , 0 , ip respectively as 
indicated in Figure 2. Under normal conditions sail aero- 
dynamic forces will create positive forces and moments on 
the boat. 

Through out the remainder of this paper the moments K, 

M, and N may be referred to interchangeably as forces or 
moments . 
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The resolution of forces into equilibrium requires that 
the summation of all the forces acting in the six directions 
be equal to zero, i.e.; 

EX. = EY. = EZ. = EK. = EM. = EN. = 0 
111111 

where " i" designates individual aerodynamic or hydrodynamic 
forces or moments acting in the prescribed direction. An 
example of a contributing force is X s , the sail driving force. 
This is one of a number of forces acting along the X q axis. 

The achievement of an equilibrium of forces in the six 
directions requires the solution of six non-linear simultaneous 
equations. To obtain a unique solution, six variables are 
required. A selection of variables for the six degrees of 
freedom problem is : 

(U b , i p, s, <J>, 9, a r ) 

The only unidentified variable is U^, which is the boat's 
velocity in the positive X q direction. Each variable contri- 
butes a major effect to a corresponding force. The boat's 
resistance is strongly dependent on boat speed, U^. As lee- 
way angle, ip, strongly contributes to side force, so does 
heel angle, <p , to the righting moment and rudder angle, a , to 
the yaw balance. 
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The present model excludes the solution in heave and 
pitch. This exclusion of two degrees of freedom was not due 
to their lack of importance, but the lack of an adequate 
parametric model. In reality the axis about which the boat 
heels may be pitched as severely as 15°. This pitch makes 
a measurable change in the boat’s angle of incidence when 
heeled to large angles. Further, due to the drop in dynamic 
pressure as the water passes around the hull a certain amount 
of sinkage will occur as a measure of the hydrostatic 
reaction. Besides the effect on form drag, this dynamic 
pressure and resulting sinkage also increases the boat's 
wetted surface and viscous resistance. The magnitude of 
these effects is small enough that they may be reasonably 
ignored. In the generation of mathematical models for the 
remaining forces, emperical data was used that came from 
experimental tests which included the effects of both 
pitch and heave. 

The present model fixes the pitch angle, 0, and sinkage, 
s, at zero. The variables are reduced to the following four: 

( U fa , ip, <p, a r ) 

The remaining major forces acting on a sailing boat are 
resolved into X, Y, K, and N components and are added together 
to give the following force equilibriums: 
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0 



1 





2 




3 




4 



The individual forces are described as follows: 
X = Sail driving force 



= Hull and keel induced drag 
X r = Rudder induced drag 

Y = Sail side force 
s 

Y. = Hull and keel side force 
k 

Y = Rudder side force 

r 

K g = Sail heeling moment 
K h = Hull righting moment 
= Keel heeling moment 
K r = Rudder heeling moment 



s 



X. 



h 



Hull resistance 
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= Yawing moment from sail driving force 



N = Yawing moment from sail side force 

N^. = Yawing moment from keel side force 

= Yawing moment from rudder side force. 

Each individual force is defined in later sections as a 
function of the boat's geometry, a true wind vector and the 
four variables U^, ifj , <£, and a^. The solution of equations 
1 through 4 for a given boat geometry and wind condition 
provides values for boat speed and attitude. From these 
values the hydrodynamic performance of the specific sailing 
boat may then be evaluated. 
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II. FORCE AND MOMENT MODELS 



2 . 1 Hydrodynamic Forces 

The hydrodynamic forces acting on a sailing boat may be 
conveniently divided into four areas for modeling. The four 
areas are hull forces, keel forces, rudder forces, and hydro- 
static forces. Hull forces consist of resistance from wave 
making, viscous effects, and form drag. Keel forces are 
made up of lift on the keel and hull, and drag induced by 
this lift. Rudder forces in the same way consist of lift 
and induced drag from the rudder. Finally, hydrostatic 
forces are caused by the change in position of the center 
of buoyancy in the hull. Each of these areas is presented 
individually and models are developed for their contribution 
to the total force equilibrium. 

2.1.1 Hull Resistance 

The estimation of hull resistance is a problem faced 
by naval architects in the design of every ocean going 
vehicle. The difference, in predicting a sailing boat's 
resistance, is in the determination of the effects changes in 
attitude have on these calculations. The approach taken is 
to determine a model for up-right resistance and modify this 
resistance as necessary to account for changes in the boat's 
attitude . 
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Upright Resistance 

The technique commonly used for prediction of resistance 
from tow tank data is based on Froude ' s Hypothesis. Basically, 
the hypothesis states that the total upright resistance may 
be divided into two components. The first is frictional or 
viscous resistance. Frictional resistance is dependent on 
Reynold's Number, a non-dimensional parameter defined as: 



R n - U b L/Y H 2 0 



The second component is residual resistance. The residual 
resistance is the combination of drag due to the generation 
of waves and the effects of eddy making and form drag. The 
effect of the hull generating waves is a gravitational effect, 
largely dependent on Froude Number. This non-dimensional 
coefficient is defined as: 



F = U, //gL 
n b = 

Though Froude left the separation of the elements of 
resistance in these two components, form drag and eddy making 
resistance have little dependence on either Froude or Reynolds 
Numbers over the velocity ranges encountered by sailing 
boats. For this reason it is logical to further separate 
residual resistance into wave resistance and form resistance. 
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The separate components of upright resistance are assumed 
independent of each other so that the total resistance of a 
boat at a given speed may be estimated by superposition of 
the forces. Total upright resistance is then: 

Upright resistance = frictional resistance 
+ wave resistance + form resistance. 

Each component of upright resistance is examined in detail 
and mathematical expressions are developed to quantify their 
magnitude . 

Wave Resistance 

The mathematical model for wave resistance is based on 
a generalized form of Michell's Integral. This expression 
is evaluated for a base hull form and is corrected for 
parametric deviation from this parent form. This expression 
is emperically adjusted to give accurate prediction of the 
towing tank results of Antiope (3) . 

Michell (7) postulated that wave resistance for thin 
ships depends upon the computation of pressure changes, 6 , 
due to the wave pattern. His expression for this resistance 
is 



1 



2 2 
(I +J ) 



a 2 -i 



dX 



5 
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where 




, 2 /r 2 

-A gz/U, , . , 

b cos ( ) dx dz 




6 




sin dx dz 
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Another form of these equations was derived by Weiblum (8) . 
A general expression of his result is 



By assuming that the block coefficient, C^, varies very little 
between yacht forms, equation 8 may be written as: 

3 

R, = —r [ d 2 f(k~, form coefficients, F ) 9 

w pg (L/ioor c n 

The number of experiments conducted to validate Michell's 
result is small. Wigley (9) compared model results to theory 
and observed the interference humps in the resistance curve 
were overpredicted by theory. Further comparison between 
vessels of constant draft and varying beam done by Wigley (10) , 
Weinblum (11) , and Mumford (12) indicate for slender ships 
wave resistance is approximately proportional to the square 
of the beam as estimated by Michell's Integral, but as the 
beam to length increases dependence of resistance on beam 
appears to be less than quadratic. It is expected that the 



R. 



B H L 

= Pg — 7 — F (sr - / form coefficients, F ) 
w jj i n 



8 
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high beam to length ratios encountered in sailing boats will 
strain the slender ship limits of Michell's Integral. 

The basic form of the wave resistance is rewritten as: 

2 

R = [ T-] n f(L/T , form coefficients, F ) 10 

w P< ? (L/100) 3 c n 

where n is between one and two. The actual value is 
determined by emperical data discussed in the following 
sections. *- 

The identification of the most important parameters 
affecting wave resistance is required. Todd (13) did a 
regression analysis of over two hundred destroyers and found 
the following coefficients had the highest correlation between 
the speed length ratios of 0.8 and 1.25. 



(V/ (L/100) , C , LCB , L/T , F ) 

P U II 



For geometrically similar boats we know directly from 

Froude's Hypothesis that the wave resistance scales is the 

cube of the characteristic length. The general form of 

Michell's Integral, equation 10, is used for the dependence 

on displacement to length ratio. By assuming base value for 

prismatic coefficient, C^; longitudinal position of the 

center of buoyancy, LCB ; and length to draft ratio, L/T ; 

o c o 

an expression for the wave resistance of a parent form may 
be written: 
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A 



11 



R = f . (U. //L) — [ =•] 

^ 1 b pg (L/100) 3 



n 



at C , LCB , and L/T 
po o co 



To predict the wave resistance of a sailing boat, the 
resistance of the parent boat is adjusted by correction 
factors based on the boat's geometric deviation from the 
parent. The wave resistance equation is 



R r = R (1 + C correction factor + LCB correction 
w wp p 

factor + L/T c correction factor) 12 

Each component of this equation is defined and discussed 
in the following sections. 

Characteristic Length - L 

The characteristic hydrodynamic length or "effective 
length" of a sailing boat is difficult to quantify. This 
has been a difficulty in the maintenance of all offshore rating 
rules as demonstrated by the extremes to which designers have 
gone to make rated length look shorter than the boat's actual 
hydrodynamic length. For the purposes of this model, L is 
determined to be the load water line, LWL. This creates 
little difficulty when examining the comparitive performance 
of similar forms, but subjective adjustment would be required 
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when comparing plumb bow, transom stern boats to long ended 
boats with shallow buttocks and entrance profiles. To some 
extent the latter might increase in length with Froude Number 
as the boat will be supported more by her overhangs. 

Parent Froude Number Dependence - f 1 (U^/ZL) 

The speed length ratio, U^/Zl, was adopted instead of 
the dimensionally correct variable F^ = U^/ZgX due to the 
common usage and the comfortable feel most yacht designers 
have with speed length ratio. The function, f^OJ^/ZL) adjusts 
the wave resistance equation to accurately predict the 
resistance of the yacht Antiope . The function was determined 
by inputting Antiope ' s geometry into the model and selecting 
f ^ (U^/ Z l) to insure proper correlation with the full scale 
towing test results reported in (3) and (14) . Values of the 
function were established at each data point given in (14) 
and a spline cubic curve was fitted to give the curve in 
Figure 3 for n = 2.0. 

3 n 

Displacement to Length Factor - [A/ (L/100) ] 

The high beam to length ratios, B/L, encountered in 
sailing boats requires that the exponent, n, be determined 
from emperical data. Unfortunately, the data base from tests 
of actual sailing boats is small. The scale effects pointed 
out in reference (15) precludes the use of results from towing 
tests of small scale yacht forms. 
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Parent Froude Number Dependence 
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Examination of Taylor Standard Series (16) at prismatic 
coefficients typical in yacht design indicate that a value 
of n between 1.8 and 2.0 appears reasonable. Though the 
prismatic coefficients, length to draft and displacement to 
length ratios in reference (16) are typical of sailing boats, 
the block coefficient, C^, is high and not characteristic of 
yacht forms. 

Figure 4 is a plot of two curves predicting the wave 
resistance of a Standfast 43 sailboat by the mathematical 
model with n equal one and two. Also plotted are the data 
points from the towing tests of a one sixth scale model 
conducted by M.I.T. and Delft University. Figure 5 compares 
the full scale test (15) of the one ton yacht Bullet to t\ra 
theoretical curves. The form drag data from both yachts was 
taken as the residuary drag at a speed length ratio equal to 
0.4 and subtracted from the residuary resistance versus speed 
curve to obtain the wave drag. 

Two is maintained as the value of n in the computerized 
model appended to this paper. Though comparisons made in 
Figures 4 and 5 indicate that n = 2 may slightly overestimate 
the effect displacement to length ratio has on wave making 
resistance, the limited sample size available for correlation 
does not justify adjustment from the theoretical value. 
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Figure 4 - Stancj-fast 43 Wave Resistance 
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Figure 5 ~ Bullet Wave Resistance 



C Correction Factor 

_JE 

The prismatic correction factor corrects the magnitude 
of the sailing boat's wave resistance for deviations away 
from the base prismatic coefficient, C p Q • The t> ase prismatic 
coefficient, C p 0 ' was selected to be the optimum prismatic 
coefficient, for each speed length ratio. Taylor's 

Standard Series (16) data are the only tests available of 
large models where prismatic coefficient was systematically 
varied. Examination of these data reveals a definite optimum 
prismatic coefficient, C p Q p t / for each tabulated speed length 
ratio. Figure 6 is a spline cubic fit of this optimum 
prismatic taken at the displacement to length ratio equal to 
250 and plotted versus speed length ratio. This optimum 
prismatic coefficient, C p Q p t ' as a ' function of the speed 
length ratio is selected as the base prismatic coefficient, 
Cp Q , of the parent equation. This is done so that deviations 
from this optimum will always result in an increase in wave 
making resistance. 

From Taylor's Standard Series, reference (16), the 
percentage increase in resistance for deviation from the 
optimum prismatic coefficient was tabulated for each speed 
length ratio. From this emperical data the expression for 

the C correction factor is generated. 

P 



Cp correction factor = f 2 (U^/ZL)] 100 ( C p 0 pt -C p^ 



1.7 



13 
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Figure 7 is a plot of the spline cubic fit of the correction 



factor's dependence on speed length ratio, f 2 (U^/ZL) . A 
reasonable definition for the sailing boat's prismatic 
coefficient, C^, must be established to insure quantitative 
agreement between the Taylor Series, reference (16) , forms and 
the yacht forms. A method for weighting the volume and 
sectional area of the keel in the calculation of the longi- 
tudinal prismatic coefficient is developed. Equations 6 and 
7 of Michell's Integral indicate that displacement should be 
weighted exponentially in the vertical direction by 



e -x 2 gz/u^. 

There are two common practices in calculating the of 

a sailing boat. The first is to include the total volume 

and sectional area of the keel into the calculation. The 

inclusion of all the deeply placed sectional area of the keel 

would overestimate its effect on the wave making resistance 

when compared to the keel's actual effect on the wave train. 

On the other hand, the second method is to treat the keel 

as an appendage and calculate only a fairbody or canoe 

prismatic. This ignores the keel's displacement effect 

altogether. Typical values of C for the same yacht calculated 

P 

both ways are 0.52 and 0.58 with and without the keel 

included respectively. If the Taylor Series, reference (16), 

optimum is to be used, the corresponding prismatic coefficient 

would be somewhere in between these two values. 
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Dependence on Froude Number , 
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One would expect the volume of the keel to have its 

largest effect when the sailing boat is sitting on a wave 

of approximately the same length as the waterline. At this 

point the wave trough would bring the keel the closest to 

the free surface. An estimation of keel's proper effect on 

prismatic coefficient, C , may be made by weighting the keel's 

P o 
- 1 /t? z 

area and volume by e nz where F is the vertical Froude 

2 nz 

number evaluated at speed length ratio equal to one. The 
longitudinal prismatic coefficient, C^, is defined as: 



c p = 



V + 
c 



v k 



(A + A* ) 
m k 



14 



where 



v» = V e -32. 23/(1. 689L) 2 



15 



and 



A* = A e ~32 . 2d/ (1 . 689L) 2 16 

k k 

V and V, are the displaced volume of the canoe body and 

C K 

keel respectively. A m and A^ are the maximum cross sectional 
area of the canoe body and the cross sectional area of the 
keel at that station respectively. The depth, d, is measured 
from the static waterline to the center of the keel's displaced 
volume . 
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L/T c Correction Factor 

Taylor Standard Series (16) tabulates data on ships at 
length to draft ratios equal to eight and sixteen. Taylor 
Series (16) asserted that linear interpolation between these 
ratios was adequate. 

The base length to draft ratio, L/T , of 16 is 

co 

selected from Taylor Series (16) data and the correction 
factor is 

L/T c correction factor = -f^ (U b //L) (16-L/T c ) 17 

where f ^ (U^/ZL) is a spline cubic fit of the speed length 
effect as shown by Taylor (16). f (U^/Zl) is plotted in 
Figure 8. 

The major differences between the ships tested for 

Taylor's Standard Series (16) and yacht forms are in the 

shape of the sections. The Taylor form typically has a 

maximum section coefficient, C , of 0.98 while sailing boats 

m 

range between 0.60 and 0.70. As a result, the vertical 

distribution of volume between the two forms differ to a 

large extent. It is not completely accurate to assume 

that vertical shifts in the volumetric distribution of the 

two forms will have a similar effect on wave making resistance. 

The model utilizes the best data available to date, but the 

need exists for a large scale yacht series examining just 

such a systematic parameter variation. 
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Figure 8 

L/ T c Correction Factor 
Dependence on Froude Number 
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LCB Correction Factor 



In the same manner as the C correction factor, the LCB 

P 

correction factor is based on the actual LCB and an optimum 
LCB position. A summation of recommended optimum LCB positions 
as a function of block coefficient and prismatic coefficient 
was presented by Todd (17) . There is reasonable agreement 
among the recommended curves with less than 1.5% variation 
between extreme curves in the area of interest. Comparing 
the LCB curves with common practice in yacht design a curve 
attributed to Professor Troost in (17) was adopted. In 
modified form it is 



LCB . = 62.5 - 17.5C 18 

opt p 

where LCB ^ is measured in percent L from the forward 
perpendicular . 

Examination of the major series in references (17) , (19) , 

and (20) which consider LCB as a variable, none adequately 
cover the region in which sailing boats fall. In all cases, 
these series cover block coefficients higher than typical 
yachts. Certain trends may be determined though: 

1. As with the previous factors, there is a peak 
effect at speed length ratio equal to one. 

2. As block coefficient decreases, the penalty for 
non-optimum center of buoyancy position also 
decreases . 

The above series, therefore, place an upper bound on the 
expected penalty. 
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